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NASA AEROTHERMAL MODELING 


OBJECTIVE: ASSESS CURRENT COMBUSTOR MOOELS 

AND IDENTIFY MODEL DEFICIENCIES 

APPROACH 

TASK I MODEL DEFINITION 

DATA BASE GENERATION 
BENCHMARK TEST CASES 

TASK II MODEL EXECUTION 

MODEL ASSESSMENT 

PROGRAM PLAN FOR MODEL IMPROVEMENT 



MODULAR APPROACH FORMS THE BASIS OF 
COMBUSTION ANALYSIS AT GTEC 


ANNULUS FLOW MODEL 

• PRESSURE DROP 

• AIRFLOW DISTRIBUTION 

• ROUNDARY CONDITIONS 


COMBUSTOR FLOW MOOEL 

• INTERNAL FLOW FIELD 

• TEMPERATURE i F/A 
DISTRIBUTION 




vz-i-s. 


TRANSITION MIXING 

• JET MIXING 

• BURNER EXIT 
TEMPERATURE QUALITY 


NEARWALL MOOEL 

• LINER CONVECTIVE AND 
RADIATIVE FLUXES 

• LINER WALL TEMPERATURE 
DISTRIBUTION 
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EMPIRICAL/ANALYTICAL DESIGN APPROACH 
HAS BEEN SUCCESSFUL IN SEVERAL 
COMBUSTOR DESIGNS 
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SUBMODELS ARE USED TO DESCRIBE 
THE FOLLOWING SEVEN COMBUSTION PROCESSES 

• TURBULENCE 

• SCALAR TRANSPORT 

• CHEMICAL KINETICS 

• TURBULENCE/CHEMISTRY INTERACTION 
(GASEOUS COMBUSTION) 

• SPRAY EVAPORATION/COMBUSTION 

• SOOT FORMATION AND OXIDATION 

• RADIATION _ 
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k-e TURBULENCE MODEL IS WIDELY USED FOR 

RECIRCULATING FLOWS 


• TRANSPORT EQUATION 

(C — D)0 = S 0 

G k = Gk (MEAN VELOCITY GRADIENTS) 
S k = G k -p< 

S e = (Cj G k -C 2 P«)j^ 
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• WALL FUNCTION ARE USED FOR NEM-WALL REGIONS 
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LOW-REYNOLDS NUMBER CORRECTION IS NEEDED 
FOR ACCURATE NEAR-WALL Z, k AND c PREDICTION 

SOURCE TERMS IN k — EQUATION = S k — 2p 

SOURCE TERMS IN t EQUATION = S c — 2 p e -° 5 y* g 

y 7 


S £ - (C-| G|< — C 2 f 2 P c) c/k 


Peff = P + Pt fp 
f 2 = 1.0 — 0.22 exp (-R t /6) 2 
fp = 1.0 — exp (— 0.011C y + ) 
y + = Pk 1 ^ C D % y/p 
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EXTRA STRAIN CAUSED BY 
STREAMLINE CURVATURE AND SWIRL 
CAN BE PARTIALLY ACCOUNTED BY 
RICHARDSON NUMBER CORRECTION 

C 2 = 1.92 exp ( - Ri V 0 ~ a c Ri c ^ <* = 0.2 
k 2 V R a 

Ri c ■ 72 m Wr) 
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ALGEBRAIC REYNOLDS STRESS MODEL 
CORRECTLY PREDICTS INDIVIDUAL COMPONENTS 

IN SIMPLE FLOWS 

Reynolds Stress Transport Equation 

^Dt — u i u j = ~ c ij + 

Assume 
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ALGEBRAIC REYNOLDS STRESS MODEL 
CORRECTLY PREDICTS INDIVIDUAL COMPONENTS 

IN SIMPLE FLOWS (CONTD) 
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C’l. ff. 0 and y are Empirical Constants 
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REYNOLDS STRESS TRANSPORT MODEL 
SOLVES SEVEN PARTIAL DIFFERENTIAL EQUATIONS 


— Diff) UjUj = S, 
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REYNOLDS STRESS TRANSPORT MODEL 
SOLVES SEVEN PARTIAL DIFFERENTIAL EQUATIONS (CONTD) 
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SCALAR TRANSPORT PROCESSES MUST BE 
MODELED PROPERLY TO ACCURATELY PREDICT 

REACTION RATES 


Assume Constant a f 

v (pVf - M vf) = S( 

Scalar Transport Equations: 

k-t Model Uses Gradient Diffusion Assumption 
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THE ALGEBRAIC SCALAR TRANSPORT MODEL 
PREDICTS COUNTER-GRADIENT DIFFUSION OF 

SCALAR IN SIMPLE FLOWS 

Transport Equations for u.fl' ; 

Div(p7 u.0') - Diff (p u.0 ? ) = PP. rt - Pc.. + p W/ . 

J J Id iff v )6 


Convection Diffusion Production | Redistribution 

Dissipation 

Assumption: 
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THE ALGEBRAIC SCALAR TRANSPORT MODEL 
PREDICTS COUNTER-GRADIENT DIFFUSION OF 
SCALAR IN SIMPLE FLOWS (CONTD) 

The transport equation for scalar fluctuations, 9 , is 

-TZ 

div (Pv# ) . Diff(P# ) = PP e - peg 

, ' * " ~~~* 

Convection Diffusion Production Dissipation 

Assume: 


— ~ |2 ' L 

div (P7d) - Diff (PQ* ) « C e (P-f) 
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THE ALGEBRAIC SCALAR TRANSPORT MODEL 
PREDICTS COUNTER-GRADIENT DIFFUSION OF 
SCALAR IN SIMPLE FLOWS (CONTD) 



Recommended Values For Model Constants: 

C-j = 1.6 C2 = 0.5, Oq - 1.25, Cff = 0.8, a-j = 0.8 

0 0 
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DATA BASE FOR IDEAL ELEMENT TESTS 
COMPILED FROM A LITERATURE SURVEY 


• DATA BASE COMPILED FOR THE FOLLOWING SUBMODELS: 

- TURBULENCE MODELING 

- GASEOUS COMBUSTION 

- SPRAY EVAPORATION AND COMBUSTION 

- SOOT FORMATION AND OXIDATION 
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DATA BASE FOR IDEAL ELEMENT TESTS 
COMPILED FROM A LITERATURE SURVEY (CONTD) 


• FOR ALL CASES, THE DATA BASE IS ORGANIZED 
IN INCREASING ORDER OF COMPLEXITY OF THE FLOW. 

FOR TURBULENCE MODELING THE CATEGORIES ARE 

- SIMPLE FLOWS (BOUNDARY LAYERS, MIXING 
LAYERS, ETC.) 

- STREAMLINE CURVATURE EFFECTS (CURVED DUCTS, 
CURVED BOUNDARY LAYERS, ETC.) 

- RECIRCULATING FLOWS (NONSWIRUNG) (BOTH 
UNCONFINEO AND CONFINED) 

- SWIRLING FLOWS (WITH AND WITHOUT RECIRCULATION) 

■ SCALAR TRANSPORT _ 
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DATA BASE FOR IDEAL ELEMENT TESTS 
COMPILED FROM A LITERATURE SURVEY (CONTD) 


• IN GASEOUS COMBUSTION, THE DATA BASE 
IS CATEGORIZED INTO 

- LAMINAR PREMIXED FLAMES 

- LAMINAR DIFFUSION FLAMES 

■ TURBULENT PREMIXED FLAMES 

- TURBULENT DIFFUSION FLAMES 

• DATA BASE FROM GARRETT GAS TURBINE 
COMBUSTORS IS INCLUDED 
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A NUMBER OF SIMPLE FLOWS 
HAVE BEEN ANALYZED 


WATTS AND BRUNDRETT 
FLOW OVER A FLAT PLATE 


EMERY ANO 6ESSNER 
2 0 CHANNEL FLOW 


EL TELBANV ANO REYNOLDS 
COUETTE FLOW 




OARIIN AND JONES 
DEVELOPING PIPE FLOW 


LAIIFER PIPE FLOW 


SAIT ANO PEERLESS 
TWO STREAM MIXING LAYER 
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A NUMBER OF SIMPLE FLOWS 
HAVE BEEN ANALYZED (CONTD) 


CHAMPAGNE AND WYGNANSKI 
MIXING OF TWO COAXIAL JETS IN AMBIENT AIR 


WYGNANSKI AND FIEDLER 
AXISYMMETRIC FREE JET 


CHANNAY ET AL 

FLOW OVER A HEATED FLAT PLATE 


Typical Axial Maan 
Velocity Profile In 
Developing Region 

Velocity Profile In 
Similarity Region 



Similarity _ 

Region 
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A NUMBER OF SIMPLE FLOWS 
HAVE BEEN ANALYZED (CONTO) 


NAUTMAN ET AL 
PLUG FLOW REACTOR 


T 


D * 10cm 


PREMIXED 

REACTANTS 


X 


PRODUCTS 


/P77777777777^777777777 

ADIABATIC WALLS 
— 120cm 


SNIPMAR ARO CO-WORKERS 
CONFIHEO STOICHIOMETRIC FLAME STARIU2ED OH A 
CTUHORICAL FLAMENOLOER IN A RECTANGULAR DUCT 


MITCHELL ET AL 
CONCENTRIC FUEL AND AIR JETS 
CONTAINED IN A CYLINDRICAL COMBUSTOR 



plow 

PATTERN 


ADIABATIC WALL 
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NASSAN Am LOCKWOOD 
FREE METHANE TUROULENT JET FLAME 
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A NUMBER OF COMPLEX NONSWIRLING FLOWS 

HAVE BEEN INVESTIGATED 


SNfVA-PRASAO AMD RAMA-PRIYAN 
FWW IN A CURVED CHANNEL 


RUNE. AMO JOHNSTON (19701 ANO EATON AM JOMSTON 
FLOW OVER A OACKWARO FACING PLANE STEP 


measurement 

LOCATIONS 



_ inlet 
straight 

SECTION 


2.54M 
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PNATARAPNOUK AND 106AN 
FLOW OVER A RIM IN A PIPE 



WOON ANO RUDIN6ER 
8U00EN PIPE-EXPANSION 
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A NUMBER OF COMPLEX NONSWIRLING FLOWS 
HAVE BEEN INVESTIGATED (CONTD) 


WED6E4MAPEO FLAMENOLDER 


API CMWSTIM TUNNEL 


CALCULATION 

DOMAIN 


air ♦ r je l ; 



COMPUTATION 

REGION 




A NUMBER OF COMPLEX NONSWIRLING FLOWS 
HAVE BEEN INVESTIGATED (CONTO) 


SCHEFFER AW SAWYER LEWIS AW SHOOT 

OPPOSED JET COHRUSTOR AX (SYMMETRIC COMHUSTOR FIT! COAXIAL FUEL AW A« JETS 
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A NUMBER OF 
HAVE BEEN 



6TEC 

NONREACTING SWIRLING COMBUSTOR 

P3-.9M ATM T-j-723°R W A 3 * /W7 W* 8 " -874 



SWIRLING FLOWS 
INVESTIGATED 

GTEC 

SWIRLING FLOW INVESTIGATION 




A NUMBER OF SWIRLING FLOWS 
HAVE BEEN INVESTIGATED (CONTD) 


BRUM AND SAMUELSEN 
SWIRL COMBUSTOR WITH COOLING AIR 



JANJUA ET AL 

SWIRLING FLOW IN A PIPE EXPANSION 
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A NUMBER OF SWIRLING FLOWS 
HAVE BEEN INVESTIGATED (CONTD) 


ALTGELO. JONES. AND WILHELM! 
CONFINED SWIRL-OR1VEN FLOW 




EL-BANNAWY ANO WHITELAW 

CYLINDRICAL COMBUSTOR WITH ROTATING CUP 
ATOMIZER AND AIR INTRODUCED THROUGH A SWIRLER 
SURROUNDING THE ATOMIZER 
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3-D DILUTION JET MIXING CALCULATIONS 
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EFFECTS OF LOW REYNOLDS NUMBER CORRECTION 
ON ASM PREDICTIONS FOR FLAT PLATE 
BOUNDARY LAYER: X = 1.8735 M 
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Y, METERS Y, METERS 

•.04 0.00 «.«4 - 0.04 0.00 0.04 


ASM PREDICTIONS FOR TWO-STREAM MIXING LAYER; 

Uf/U| = 0.43 



U-VEL, M/SEC 



U-VEl. M/SEC 


-X-IT.S CM 
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ORIGINAL PAGE ESI 
OF POOR QUALITY 


ASM PREDICTIONS FOR 
TWO-STREAM MIXING LAYER, 
UE/UI = 0.43 









PREDICTIONS FROM ASM WITH 
STREAMLINE CURVATURE CORRECTION 
ALONG CONVEX (INNER) WALL 


X-.2SIS METER* 


X-.47S0 METERS 


ORIGtNM- Wj** 
of POOR Q uWJTT 


UV, M2/S2 

X* .7010 METERS 


UV, M2/S2 

X-.0240 METERS 


UV, M2/S2 

I- 1.14* METERS 


UV, M2/S2 

1.00 

UV, y2/S2 

2.00 






OKiuHVAL W 

OF POOH QUALfTY 


ASM PREDICTIONS AND DATA FOR 
(OUTER) CONCAVE WALL 
ROUNDARY LAYER 


X-.t9l« MfTSM 


X-.4TIO NITERS 


UV, M2/S2 

X-.TOIO NITERS 


UV, M2/S2 

X-.I240 NITERS 


UV, M2/S2 

UV, y2/S2 

X-1.140 NCTIRS 
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1.00 2.00 
UV, U2/S2 





PREDICTIONS BY ASM WITH STREAMLINE 
CURVATURE CORRECTIONS FOR (OUTER) 
CONCAVE WALL BOUNDARY LAYER 
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2.00 
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COMPARISON BETWEEN DATA AND PREDICTIONS 
BY ASM WITH STREAMLINE CURVATURE 
CORRECTIONS FOR FLOW BEHIND 3.81 CM STEP 
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STANDARD K-e MODEL PREDICTIONS FOR PREMIXED 
PROPANE/AIR GAS FLOW BEHIND 2.54 CM STEP; 

0 = 013 


2-STEP KINETIC SCHEME 
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4-STEP KINETIC SCHEME 
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COMPARISON BETWEEN DATA AND 
K* MODEL PREDICTIONS FOR FREE 
SWIRLING JET; SWIRL NUMBER 

= 025 






AXIAL VELOCITY < M/S I 


Ml 
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COMPARISON BETWEEN DATA AND 
K-t MODEL PREDICTIONS FOR FREE 
SWIRLING JET; SWIRL NUMBER 

= 025 


N 





TANG. VELOCITY ( M/S ) TANG. VELOCITY (M/S) 
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COMPARISON BETWEEN DATA AND PREDICTIONS 
USING 20,000 NODES FOR J = 25.32, 

S/D = 2.00, H/D = 8X10 



X/H s 0.80 X/DJ *•!• 



X/M * 0.78 X/OJ *7.75 



X/H « 14D X/DJ «!0*J9 



C TNNllM l/C TNtUM«JI 
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COMPARISON BETWEEN DATA AND PREDICTIONS 
USING 6000 NODES FOR J = 25.32, 

S/D = 2.00, H/D = 8.00 


PREDICTIONS MEASUREMENTS PREDICTIONS MEASUREMENTS PREDICTIONS MEASUREMENTS 



X/H * 0.50 X/DJ dSAB X/H x 0.75 X/OJ *7.75 X/H s 1.00 X/OJ =10.33 








k-£ MODEL AND ITS MODIFICATIONS 


• GIVES GOOD CORRELATION OF MEAN VALUES FOR 

■ SIMPLE FLOWS 

■ FAR-FIELO REGIMES OF NONSWIRLING/SWIRLING FLOWS INVOLVING REGIONS 
OF RECIRCULATION 

■ NONRECIRCULATING SWIRLING FLOWS 

■ OUTER REGIONS OF STRONG SWIRLING FLOWS 

• GIVES REASONABLE CORRELATION OF MEAN VALUES FOR 

> NONSWIRLING RECIRCULATING FLOWS EXCLUDING VICINITY OF REATTACHMENT POINT 

■ CONFINED DISK FLOW WITH A CENTRAL JET 

■ SHEAR LAYER OF STRONG SWIRLING FLOWS 

> CONFINED SWIRLER WITH HUB AND SHROUO EXPANSIONS 


• PREDICTS TRENDS IN MEAN VALUES FOR 
> RECIRCULATION ZONE OF SWIRLING FLOW 

■ CONFINED SWIRLER WITHOUT OUTER EXPANSION 

• REQUIRES DIFFERENT MODIFICATIONS FOR CONVEX AND CONCAVE 
WALL SHEAR LAYERS 

^ 63-0205 29 A 




ALGEBRAIC STRESS MODEL: 
REYNOLDS-STRESS PREDICTIONS 

GIVES GOOD CORRELATION FOR 

- SIMPLE FLOWS 

> NORMAL STRESSES IN NONSWIRLING RECIRCULATING FLOWS 

GIVES REASONABLE CORRELATION FOR 

> SHEAR STRESSES IN NONSWIRLING RECIRCULATING FLOWS 

- NORMAL STRESSES IN SWIRLING FLOWS 

PREDICTS TRENDS IN 

- SHEAR STRESSES IN SWIRLING FLOWS 

DIFFERENT MODIFICATIONS REQUIRED FOR CONVEX AND CONCAVE 
WALL SHE1R LAYERS 

MEAN FLOW FIELD PREDICTIONS ARE SIMILAR TO K-e RESULTS 



SCALAR TRANSPORT MODELS 
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K-e MODEL WITH PRANDTL/SCHMIDT NUMBER - GOOD 
GRADIENT DIFFUSION APPROXIMATION IS VALID 


• ALGEBRAIC SCALAR TRANSPORT MODEL 
■ PROMISING APPROACH 
" MORE VALIDATION WORK IS NEEDED 


WHERE 





TURBULENCE/CHEMISTRY INTERACTION 


BOTH TWO-STEP AND FOUR-STEP SHOW PROMISE 


MODIFIED EDDY BREAKUP SHOULD CONTINUE TO DEVELOP 
BECAUSE IT CAN BE EASILY EXTENDED TO MULTISTEP SCHEMES 


BILGER’S TWO-REACTION ZONE MODEL GIVES GOOD RESULTS 
FOR JET FLAMES, REQUIRES MORE WORK 



DILUTION JET MIXING 


SLIGHTLY UNDERPREDICTS JET PENETRATION 
AT LOW TO MODERATE J VALUES 

CENTERLINE TEMPERATURES PREDICTED WELL 

TRANSVERSE MIXING PREDICTIONS SLOWER THAN DATA 

EFFECT OF S/D. H/D. J ON MIXING PREDICTED CORRECTLY 

GOOD COMPARISON WITH DATA FOR JET INJECTION FROM 
> ONE WALL 

■ BOTH WALLS - INLINE ORIFICES 
- BOTH WALLS - STAGGERED ORIFICES 



CASES WERE DIVIDED INTO FOUR FLOW CATEGORIES 

THE MAIN OBJECTIVE OF THE NASA-SPONSORED AEROTHERMAL MODELING 
PROGRAM. PHASE I WAS TO ASSESS CURRENT AEROTHERMAL SUBMODELS USED 
IN THE GARRETT TURBINE ENGINE COMPANY (GTEC) ANALYTICAL COMBUSTOR 
MODELS. 

A NUMBER OF “BENCHMARK" QUALITY TEST CASES WERE SELECTED AFTER AN 
EXTENSIVE LITERATURE SURVEY. THE SELECTED TEST CASES. BOTH NONREACTING 
AND REACTING FLOWS. WERE BROADLY DIVIDED INTO THE FOLLOWING 
CATEGORIES: 

• SIMPLE FLOWS 

• COMPLEX NONSWIRLING FLOWS 

• SWIRLING FLOWS 

• DILUTION JET MIXING IN CONFINED CROSSFLOWS 
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TURBULENCE AND CHEMISTRY MODELS 

WERE ASSESSED 


THESE TEST CASES WERE USED TO ASSESS THE FOLLOWING SUBMODELS 
SEPARATELY AND JOINTLY FOR VARIOUS COMBUSTION PROCESSES: 

• k-e MODEL OF TURBULENCE ANO ALGEBRAIC STRESS MODEL, 

WITH AND WITHOUT VARIOUS CORRECTIONS INCLUDING LOW 
REYNOLDS NUMBER AND RICHARDSON NUMBER CORRECTIONS 

• SCALAR TRANSPORT MODELS 

• MULTISTEP KINETIC SCHEMES 

• TURBULENCE/CHEMISTRY INTERACTION 

• SPRAY COMBUSTION 







ADVANCED NUMERICAL SCHEME IS REQUIRED 


THE FOLLOWING GENERAL CONCLUSIONS WERE DERIVED FROM PHASE I WORK 


• AN ACCURATE NUMERICAL SCHEME SHOULD BE DtVELOPEO 
TO MINIMIZE NUMERICAL DIFFUSION IN THE COMPUTATIONS 
OF RECIRCULATING FLOWS 


• BENCHMARK QUALITY DATA SHOULD BE GENERATED UNDER 
WELL-DEFINED ENVIRONMENTS FOR VALIDATING THE VARIOUS 
SUBMODELS USED IN GAS TURBINE COMBUSTION ANALYSIS 


Ln 

o 





MORE MODEL DEVELOPMENT IS NEEOEO 



